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The Schrodinger integral-equation approach for calculating the classical first-passage time (C-fpt) 
probability density is extended to the case of quantum first-passage time (Q-fpt). Using this ex- 
tension, we have calculated analytically the Q-fpt probability density for a class of few-site/state 
tight-binding Hamiltonian systems, e.g., a qubit, as well as for an infinite ID lattice. The defining 
feature of such a quantum system is that the passage across the boundary between a subspace (w) 
and its complement (u>) is through a unique pair of door-way sites such that the first departure from 
(arrival at) uj corresponds to the first arrival at (departure from) uj. The Q-fpt probability density so 
derived remains positive over the time interval in which it also normalizes to unity. These conditions 
of positivity and normalization define the physical time domain for the Q-fpt problem. This time 
domain is found to remain finite for the few-site/state Hamiltonian systems considered here, which 
is quite unlike the case for the diffusive C-fpt problems. The door-way sites and the associated 
Q-fpt probability density derived here should be relevant to the inter-biomolecular/nanostructural 
electron transport phenomena. 
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For a particle restricted to lie initially at a point v in- 
side a finite subspace u) bounded by the surface cr, the 
first-passage time (fpt) is defined as the time of its first 
crossing of the boundary a (concomitant with its first ar- 
rival at the complementary subspace u>) [THB]. Thus, for a 
classical (C) stochastic system evolving probabilistically, 
one can conveniently introduce the restricted (r) proba- 
bility Pf?(ui,t | v>, 0) that the particle remains confined 
to the subspace u> for all times t' with < t' < t. The 
associated C-fpt probability density is then necessarily 
given by P£(o>, * | u, 0) = - (^P r c (u;, t | u, 0)) 0. 
The parameters uj and a are to be chosen as appropriate 
to the specific physical realization of the problem. (It is 
to be noted here that the crossing of a boundary may, 
in general, connote crossing any sharply defined physical 
condition, e.g., of the first departure from, or the first ar- 
rival at the subspace of interest). The fpt-probability dis- 
tribution is of interest for a number of physical, chemical 



The essential logic underlying the above integral equa- 
tion is simple: The first part of r.h.s. is the restricted 
probability of never venturing beyond the subspace uj up 
to time t. The second part is the complementary prob- 
ability that takes into account the multiple excursions 



and biological rate processes involving, e.g., the classical 
probabilistic (Kramers) escape over a potential barrier 
[2], the quantum tunneling through an activation bar- 
rier and the associated delay [7HTB]. the single-electron 
quantum transfer in nanostructures [TB] and biomolecules 
[17] . and, of course, for the classic case of nuclear (3- 
decay. For a classical stochastic process, e.g., physical 
diffusion, the fpt problem is well posed, and the problem 
has been solved variously and exhaustively PUS], e.g., by 
use of a perfect absorber -a boundary condition that can 
be implemented mathematically through the well known 
Kelvin method of images [TJ [2] • 

An elegant approach to the classical fpt-probability 
density for a ID continuum is the one based on the in- 
tegral equation relating the restricted (r) and the un- 
restricted (u) probabilities (Pf? and P£ ), as proposed 
by Schrodinger @] (well before the appearance of the 
Schrodinger wave equation of quantum mechanics), writ- 
ten symbolically as 



(1) 

I 

across the boundary a taking the particle to uj at some 
time < t' < t, and then bring it back to uj within the 
remaining time (t — t'). This is an exhaustive reckoning 
of all the possible alternatives subsumed in the P^ on 
the l.h.s of eq. 0. 



P?{w,t | i/.O) = P?(uJ,t | „,0) - J* (I^W I v,Q)) P?{w,t\ BJW 
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FIG. 1. Classical system of a particle hopping between the two 
sites 1(= lj) and 2(= Q); the dotted line depicting notionally 
the boundary a, with p the hopping probability rate. 



Thus, e.g., for a classical 2-site hopping model with 
the subspaces (or rather the subsets) ui = 1 and Co = 2 
as depicted in Fig. [T] the integral equation reads 



*f(M|l>0) 



P r c (l,t | 1,0) 



J*dt' PZ(l,t\2,t') ^Pf(l,t'|l,0). 



(2) 



The unrestricted probabilities involved here are noth- 
ing but solutions of the master equation for the system 
under appropriate boundary conditions. With this, it is 
straightforward to solve the integral equation, e.g., by a 
simple Laplace-transform [5]. The time of first passage 
then refers to the time of first departure from the initial 
site 1, which is, of course, concomitant with the time of 
first arrival at site 2. The integral equation as written 
for the case of C-fpt in eq. sets the notations for the 
quantum case (Q-fpt) henceforth for a class of many-site 
systems where the boundary (cr) separating the subspace 
to and its complement Co is notionally a single point as 
depicted by the dotted vertical line in Fig[l] 

Before we turn to Q-fpt, let us note that the Q-fpt 
problem has been generally viewed as ill-posed, and has, 
therefore, remained unsolved so far to the best of au- 
thors' knowledge (see however [S]). Reasons being: (a) 
inapplicability of Kelvin's method of images inasmuch as 
quantum mechanics does not admit perfect absorber sans 
reflection UHl US] ; (b) essential difference between the 
Wiener (real) functional integral [20] of classical stochas- 
tic mechanics and the Feynman (complex) path integral 
[21) of quantum mechanics; (c) non-demolition unobserv- 
ability of quantum paths. (However, see [55] [53] for the 
possibility of weak measurements); and finally, (d) the 
question of time itself being an operator [24H27] . (How- 
ever, see [28]). 

In any case, the Q-fpt probability distribution is op- 
erationally well defined -after all, for any physical event 
that may come to a pass, there must necessarily be a time 
of its first passage past the post (i.e., the marker or the 
boundary). To us, this seems to be a robust idea. The 



question then really is how it is to be calculated quantum 
mechanically. It is this central issue that we turn to now. 

For the class of discrete (N-site) problems with only 
nearest-neighbor couplings as considered here, the sub- 
space (or rather the subset) u> has a unique single site of 
departure to its complement Co, and Co has a unique single 
site of arrival from oj.(It is apt to call such a unique pair 
of sites bridging to and Co as the door-way sites in gen- 
eral). Now, this uniqueness of connectivity eliminates 
any quantum interference of alternatives, -other than 
those already contained in the unrestricted probabilities 
P$(l,t | 1,0) and P$(l,t | 2,1/) occurring in eq. 
Thus, the Schrodinger integral equation for the classi- 
cal case remains valid even for the quantum case, but 
with the proviso that the two unrestricted probabilities 
referred to above are calculated quantum mechanically. 
The restricted probability P® is then determined implic- 
itly through the integral eq. ([2| itself. Thus the problem 
of quantum phase in connecting the restricted and the 
unrestricted probabilities at the point of return form Co 
in eq. ^ disappears. With this proviso, we now turn to 
calculating the Q-fpt probability density. 

First, we consider the two site/state quantum sys- 
tem (essentially a qubit) described by the tight-binding 
Hamiltonian 

n 2 = e 1 |l)<l| + e 2 |2>(2|-ft 7 |2)(l|-ft 7 *|l)<2|, (3) 

where |1) is taken to correspond to the subspace u>, and 
1 2) to the complement subspace Co. Here for simplicity 
of calculation, we will set the site energies e\ = = e 2 , 
and take 7, the tunneling matrix element, to be real with 
7 = 1 = h. Here, the boundary (a) is notionally a single 
point as depicted by a vertical dotted line in its classical 
analogue in Fig. [T] . For this quantum two-site system, 
the Schrodinger integral equation reduces to 



i£(M|l,0) 

= P®(l,t 1 1,0) 



jf'df i*(M|2,f) ^P?(l,t'\i,0). 



(4) 



The structure of eq. Q embodies the following physi- 
cal conditions: 

(A) . The restricted probability P,9(l,t|l,0) decreases 
monotonically from unity at t = to its first zero 
at some time T, which can even be infinite; 

(B) .It follows from (A) that the Q-fpt probabil- 
ity density Pg(l,t|0,0) = -f P?{1, t \ 1, 0) 
is positive (but not necessarily monotonic), and 
is automatically normalized to unity over the 
time interval < t < T , as J Q T -§- t P®{l,t | l,0)dt 
= P r Q (l,0 | 1,0) - P r Q (l,T I 1,0) = 1 ; and, 

(C) .It follows from (A) and (B) that the time interval 
< t < T is the physical time domain of definition 
for the Q-fpt problem. 
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FIG. 2. Quanutm 2-site system: plots of the restricted prob- 
ability P?(l, t | 1,0) (dashed line) and the first passage 
probability density P® p (l,t j 1,0) (solid line) against t, with 

T = 275' 
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FIG. 3. Quanutm 3-site system: Plots of the restricted prob- 
ability P®(u), t \ 1,0) (dashed line) and the first-passage prob- 
ability density P^ p (oj,t j 1,0) (solid line) against t, with 
T = 1.72757. 



The two unrestricted probabilities in eq. Q are ac- 
tually the results of the quantum unitary evolution as 
described by the time-dependent Schrodinger equation. 
They are obtained as follows : 

P«(l,t|l,0) = |(l| e -^ t |l)| 2 =cos 2 (i), (5) 
PQ (l,t | 2,0= Klle-^ 4 -''^)! 2 = sin 2 (i- t') . (6) 

This is, of course, not the case for the restricted proba- 
bility P®{t). But, as discussed earlier, the whole point 
of our treatment is that this restricted probability (t) 
gets determined implicitly through the Schrodinger inte- 
gral equation Q in terms of the unrestricted probabilities 
(P^p) for the class of quantum systems with the door-way 
sites as considered here. 

Substituting from eqs. ^ and ([6| into eq. Q yields, 
after some algebra involving Laplace transforms, the re- 
quired solutions as given below and plotted in Fig. [2] 



P r Q (M | 1,0) = cos(v / 2i), 
P%(l,t | 1,0) = y/2 sin(V2t). 



(7) 
(8) 



The above treatment can be readily extended to the 
other few-site cases with appropriate choices for the sub- 
spaces to and uj. Thus, for the 3-site quantum system 
{|1), 1 2), 1 3)} with the tight-binding Hamiltonian 



H 3 = -|1)<2|-|2)(1|-|2)(3|-|3)(2| 



(9) 



let us consider the case in which ui = {|1),|2)} and 
Cu = 1 3). When the system is initiated at time t = 
at site |1) of w, the following relevant unrestricted prob- 
abilities are obtained: 



pQ(u,t | 1,0) = Klle^ll)! 2 + K2|e-* M » 1 |1> 



3 — cos ( -y/2 t 



cos 



V2 



(10) 



FIG. 4. Quanutm 4-site system: Plots of the restricted prob- 
ability P®(uj, t | 1,0) (dashed line) and the first-passage prob- 
ability density P^ p (oj,t | 1,0) (solid line) against t, with 
T= 1.86574. 



P®{uj,t | 3,0 = Kl|e- <w »< t - t '>|3)| 2 + |(2|e-^ 3 ( t - t ')|3)| 2 
3 + COS (-y/2 (t - t')) . 2 {t-t n 



V2 



(11) 



In terms of these, we have the exact solutions given below 
(with the numerical coefficients rounded off) and plotted 
in Fig(3| 

P r Q (w,i | 1,0) = 1.132 cos(0.915 t) - 0.132 cos(2.676 t), 

(12) 

P?Ju,t | 1,0) = 1.036 sin(0.915 t) - 0.354 sin(2.676 t). 

. (13) 
Similarly, the results obtained for a 4-site system 

{|1), |2), |3), |4)}, with w = {|1), |2)}; Q = {|3), |4» and 

v = |1) are given below and plotted in Fig|4j 

P, Q (w,t | 1,0) = 0.954 cos(0.754 t) 

+ 0.131 cos(1.261 t) - 0.085 cos(2.973 t), 

(14) 
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P? p (u,t | 1,0) = 0.719 sin(0.754 t) 

+ 0.165 sin(1.261 t) - 0.253 sin(2. 973 t). 

(15) 

It should be noted here that for all the few-site cases 
treated above, the three physical conditions (A, B and 
C) listed previously are indeed satisfied, and the time T 
defining the domain for the Q-fpt problem is found to 
remain finite. Also, by comparing the results obtained 
for the 4-site and the 3-site problems, it is observed that 
the addition of extra distant sites in the complementary 
subspace Q does indeed modify the Q-fpt distribution. 
This is quite contrary to the case of diffusive C-fpt 
where it does not. This is essentially a manifestation of 
the so-called non-locality of quantum mechanics, orches- 
trated by the virtual excursions across the boundary (cr) . 

Finally, we consider the case of an infinite ID lattice 
with the tight-binding Hamiltonian 

oo 

-Hoo = + + l"X ri + 1 l)> (16) 

n— — oo 

where we choose lu = {0, —1, —2, —3 . . . , — oo} and U) = 
{1, 2,3,..., +oo}. The system is initiated at site |0) at 
time t = 0. In this case, the Schrodinger integral equation 



P®(Lj,t\0,0) = P?(w,t I 0,0) 

dt' pQ(u,t\i,t') ApO( Wjt | ,o), 

(17) 



reduces to give simpler expressions for the Laplace- 
transforms (£) of P® and Pf p : 



£[P?(u>,t | 0,0)] = 



c[pQ(oj,t\o,o)-pQ(cj,t 


1,0)] 


1 - sC 


P&{u,t | 1,0) 





c 



Pg(w,t|0,0) 



l-s£ 




o,oy 


1 - sC 


~P${u,t 


1,0) 



(18) 
(19) 



Now, with the help of the time-Laplace and the 
lattice-Fourier transforms, we solve for the un- 
restricted unitary evolutions of the wavefunction 
I^W) = Y^L-od a n(t)\ n )i given by the time- dependent 
Schrodinger equation 



-id n {t) = a n -i(t) + a n+1 (t) 
for the two initial conditions: 

(I) for o„(0) = S nt o giving 

— OO 

pQ( w ,t|o,o) = ^|4(t)| 2 

n=0 



(20) 



(1 + J 2 (2t)) 



(21) 




FIG. 5. Quantum infinite lattice system: Small time behav- 
iors of the restricted probability P®(oj,t | 0,0) (dashed line) 
and the first passage probability density P9 (u,t | 0,0) (solid 
line) against t for small time. 



and, (II) for a„(0) = 5 n ,i giving 

i*( W ,t | 1,0) = g K(t)\ 2 = {1 -f 2t)) . (22) 



Here the superscripts (I and II) on the amplitudes a n (t) 
denote the corresponding boundary conditions. Upon 



substituting from eqs. (21 1 and (22) back into eqs. (18) 



and (19), we obtain 



4 K(-i^) 

£ \P?{u>,t | 0,0)] = V f>^ , (23) 

L r v 1 ' } s 7r + 2K(-i|) v ; 



£ 



Pg(w,t|0,0) 



*-2K(-jg) 
7r + 2K(-i|) 



(24) 



where, K(x) is the complete elliptic integral of the 
first kind, and, we have used the identity £ [</q(ci£)] = 

2 k(-^t) nasni- 



Equations (23) and (24) cannot be written in terms 
of elementary functions. From the large s behavior of 
the Laplace transforms in these equations, however, we 
have obtained the small-time behavior of the required re- 
stricted and first-passage probabilities P® and Pj (using 
a professional software package like [3D]). These are plot- 
ted in Fig. [5] As can readily be seen from the plots, that 
in this limited small time regime, P^? remains positive 
while falling monotonically. Also, P? remains positive, 
though not monotonic. We emphasize that both of these 
trends are physically valid. 

In conclusion, we have extended the Schrodinger 
integral equation for the classical first-passage time 
(C-fpt) probability density to the case of quantum first- 
passage time (Q-fpt), and solved it exactly for a number 
of few-site model systems described by tight-binding 
Hamiltonians, including the case of the two-site/state 
(qubit) system. In these model quantum systems, for 
the first passage across the boundary separating the two 
subspaces/subsets, there is a unique site of departure 
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(arrival) on one subset - connecting to the corresponding 
unique site of arrival (departure) on the complementary 
subset - the door-way sites/states. This crucial feature 
takes care of the essential quantum phase problem 
involved here. The calculated Q-fpt probabilities satisfy 
the necessary physical conditions of positivity and nor- 
malization to unity over the time-domain of definition 
of the problem, which turns out to be finite unlike the 
case of C-fpt. The case of an infinite ID tight-binding 
Hamiltonian system is also treated. We hope that 
the present approach opens up a novel view upon this 
important field of Q-fpt problems whose solutions have 
so far eluded us. More specifically, this approach should 
be of special relevance to the quantum phenomena 
of intcr-biomolecular/nanostructural electron transfers 
[THIH?], where such door-way states are expected to exist. 
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